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Abstract 



We study quenched and dynamical pion correlators in the two-flavour 
Schwinger model (2D QED) both within topological charge sectors and 
averaged. The results show that a clean-cut definition of a pion mass is no 
longer possible if tunneling between sectors does not occur. 

0> PACS numbers: 11.15.Ha, 12.20.-m, ll.10.Kk 

6-i. 1 Introduction 

Constructing decorrelated configurations in lattice simulations can be difficult 
if large energy barriers exist between regions of the configuration space to be 
sampled. One such possibility are topological sectors for the U(l) model. While 
local updates do not work well, we can use in two dimensions a global heatbath 
update to obtain reasonably high tunneling rates. This trick obviously works 
also in the quenched case, while for full dynamical simulations global heatbath 
updates are not known yet. 

A way of dealing with such systems which is sometimes suggested |IJ is to stay in 
one fixed topological charge sector and define quantities in that way. In order to 
study this for relevant observables, we measure pion correlators for the massive 
two-flavour Schwinger model (QED in two dimensions) |2| in fixed (low) topolog- 
ical sectors. In the quenched case we compare local and global update schemes. 
For dynamical fermions we compare low and high (3 results, using both the local 
bosonic algorithm (LBA) [g] and Hybrid Monte Carlo (HMC) ||. 

*supported by DFG research grant No. WO 389/3-2; email: elser@linde.physik.hu-berlin.de 
t email: bunk@linde. physik.hu-bcrlin.de 



1 



2 Model 



The use of the Schwinger model with two degenerate flavours of Wilson fermions 
is motivated by the fact that it is rather QCD-like, but still tractable with small 
computational effort. It exhibits an axial vector anomaly, light pseudoscalar 
mesons (pions) and massive particles (a ,r]) 0. 
The lattice model is defined via the action 

S = S G (U) + S F (U,V), 

S G = P'EQ.-RbUp), 
p 

X fl 

where Sq is the standard plaquette part and Sf the Wilson fermion action. Up 
denotes the plaquette variable, U x ^ the compact link at site x in direction \i and 
\1/ the two-component fermion field. We use periodic boundary conditions for link 
variables and anti-periodic ones for fermions throughout. 

3 Algorithm 

Quenched case: local update. We use a local link update consisting of one exact 
heatbath |5| and three over-relaxation steps per trajectory. 
Quenched case: global update. Alternatively, it is possible to use a global heat- 
bath for the plaquettes. As the new configuration is constructed without recursion 
to the old one, we expect to encounter much less problems with metastabilities. 
The global update uses the fact that nearly all plaquettes are independent even 
on a finite lattice. The plaquettes have only to satisfy the constraint 

V] Up = 27m; n ETL. 
p 

We are therefore able to update LT — 1 plaquettes with a heatbath algorithm 
and then choose the last one using a Metropolis acceptance step. This plaquette 
configuration has to be translated into a valid link configuration. To achieve this, 
we utilize the freedom to choose a gauge. In our case we use a maximal tree 
prescription, setting LT — 2 links to 1. Then LT links can be recursively fixed 
from the plaquettes. The unconstrained two remaining links correspond to the 
free Polyakov loops Pt and Pl in 2 dimensions. 

We want to state that the problem of slow topological charge fluctuations could 
also be solved by explicit topological updates [0. Unfortunately, this trick is not 
applicable in the presence of dynamical fermions. 
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Full dynamical simulations. In this case we use the Hermitean version of Liischer's 
local bosonic algorithm (LBA) 0. A noisy Metropolis step is used to make the 
algorithm exact as described in [BL 0]. To have an independent check we have 
also implemented a Hybrid Monte Carlo (HMC) algorithm. 



4 Topological sectors 

The integer-valued topological charge functional 



47T 

can be represented on the lattice by 



with plaquette angle 4>p = Imln(JTp) £ (— 7r, tt) [l(|. We denote by tunnel events 



all updates resulting in a change of the topological charge and as tunnel proba- 
bility the number of tunnel events divided by the total number of updates. The 
topological susceptibility is defined via 

xtop = -tH< Q 2 > - < Q > 2 ]- 

JV p 

We demonstrate the relevance of topological sectors on a 16 x 32 lattice showing 
in Fig. the tunnel rate plotted against (3 for local and global updates. The 
exponential decrease of the tunnel probability for this local update algorithm 
gives rise to metastabilities in simulations at large j3. We therefore consider it 
worthwile to investigate observables within fixed topological sectors. 



5 Simulations 

Simulation parameters. We simulate on 8 x 20 and 16 x 40 lattices at a beta 
value of j3 — 12 generating 10000 configurations. We only show data for the 
larger 16 x 40 lattices. For the fermion part, we choose k = 0.24 where the pion 
correlation length (in the quenched case) is found to be around 3 and finite size 
effects can be expected to be small. 

Local updates. We monitor the topological charge during our simulations. Due 
to metastability, no tunneling of Q is observed in the runs with local updates. 
At (3 = 12 we therefore are able to perform a simulation in a given topological 
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Figure 1: Tunnel probability as a function of (3 

U( 1 ) pure gauge model 16x32 
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sector by using an initial configuration with this particular charge. This is done 
generating a classical homogeneous plaquette configuration of the desired charge 
and converting this to the links as described in Sect. ^| The two Polyakov loops 
Pt, Pl are chosen according to a flat random distribution. 

Global updates. In the limit of independent plaquettes the topological suscepti- 
bility can be calculated analytically. An approximation for f3 — > oo gives 

1 

Xtop ~ 4^3" 

which for (3 = 12 amounts to Xtop = 0.0022. To check our global update, we 
simulate on a 16 x 32 lattice obtaining xtop = 0.0021(1). This shows that the 
global update works well without metastabilities in the topological charge. 
Observables. We measure the pion correlator using the operator tp^^Tip. We apply 
a point source and wall sink prescription to increase statistics ||. The effective 
masses are calculated from two adjacent time slices using the hyperbolical cosine 
formula. The error analysis of the effective pion mass takes covariances and 
auto-correlations of the correlators into account. 

5.1 Quenched results 

The results of the quenched runs are shown in Fig. ^|. Obviously, for Q = 0,4 
we are not able to find a plateau in the effective mass plots. We find a valley-like 
structure in the mass for the low Q cases. For high Q this turns into a hill-like 
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Figure 2: Quenched effective pion mass as a function of time for (3 = 12.0, 
k = 0.24 
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Figure 3: Dynamical effective pion mass as a function of time for (3 = 12.0, 
k = 0.24 
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structure with the peak situated at half of the temporal lattice extent. To show 
that in the intermediate region the valley and hill structures can approximately 
cancel and suggest a fake plateau, we include the Q — 1 plot. 
In the quenched case, we are able to compare to the results using a global update 
scheme. For the global update we find a tunnel probability of P = 0.76 and 
therefore do not expect any influence of topology. The effective mass is shown in 
Fig. [|d. A plateau is clearly more reasonable than in the fixed Q cases. 

5.2 Dynamical fermion results 

Effective masses from full dynamical simulations are shown in Fig. |^. We do not 
expect quantitatively the same results as in the quenched case. But the behaviour 
is qualitatively similar. 

To establish that our dynamical results are not influenced by the chosen param- 
eters of the local bosonic algorithm, we repeat the calculation with topological 
charge Q = 4 using a standard HMC algorithm. This is also shown in Fig. |^. 
The results agree nicely within errorbars. 

From these results, we conclude that we need to average over the topological 
sectors to obtain a plateau in the effective mass. 

6 Projections to topological sectors 

To gain further insight, we now use a slightly different approach. In principle, we 
could also restrict ourselves to definite topological sectors by selecting measure- 
ments with fixed topological charge from a simulation, i.e. effectively simulating 
the path integral given by 

Z[q] = J D\U]D[*]D[*]8 Q jr 8 . 

To this end we need simulations with a reasonably high fluctuation. Such sim- 
ulations can be done e.g. with dynamical simulations at low (3 or quenched 
simulations using global updates. 

Full dynamical case. We work at low (3=1 with a slightly smaller k = 0.22. 
Effective masses are depicted in Fig. |[ We can detect no discrepancy between 
masses calculated in different topological sectors. This result was also reported 
by a group working with staggered fermions, which concluded that topological 
sectors to not matter for the Schwinger model [|IJ. 

Quenched case. Here we exploit the opportunity to use the same parameters (3 = 
12, k = 0.24 as in Sect. [|. The results are shown in Fig. |[ At this (3, we do 
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Figure 4: Dynamical projected effective pion mass as a function of time for (5 
1.0, k = 0.22 
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Figure 5: Quenched projected effective pion mass as a function of time for (5 
12.0, k = 0.24 
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not find agreement. Rather the effective masses are nearly the same as in the 
simulations without tunneling presented in Fig. || We would like to point out 
that they do not agree within errorbars. On the other hand, we remark that the 
statistical sample was very much smaller for the projected data due to the fact 
that only a part of the generated configurations is projected into the appropriate 
sectors. 

The striking difference between low and high (3 results makes a sound understand- 
ing of this behaviour highly desirable. It is evidently not just the averaging over 
the topological sectors (as found in Sect. §) which is lacking in high (3 simulations, 
but there seems to be some more subtle dynamical effect involved. 

7 External gauge configurations 

In order to further investigate how much averaging is necessary we plot in Fig. |] 
effective masses for external configurations with fixed topological charge Q = 
and Q = 4. These are generated from homogeneous plaquette configurations in 
the way described in Sect. 0. For the fermions we use k = 0.24. 
We stress that we use random Polyakov loops P?, Pl in both cases, so that one 
should not expect free fermion behaviour in the Q = case. For Q = 4 a 
completely irregular behaviour is observed. It is thus not possible to measure a 
meaningful pion correlator from one (even very smooth) configuration alone. The 
translation invariance is manifestly broken even for that smooth configuration. 
The result of averaging over 10 values of Pt and Pl is shown in Fig. [F[ We clearly 
regain the qualitatively expected regular valley and hill structure observed in Fig. 
| in both the Q = and Q = 4 case. 

8 Conclusion 

Our results for effective pion masses for simulations at fixed topological charge at 
high (3 clearly show that definition of a pion mass from a plateau is not possible for 
both quenched and dynamical simulations even for vanishing topological charge. 
A comparison with quenched global updates exhibiting no metastabilities demon- 
strates that a plateau can be found in the correct path integral sample. Further- 
more, cross-checks against HMC indicate that this problem is not an artifact 
stemming from the fermion algorithm. 

In Sect. H we studied effective pion masses calculated from projections to fixed 
topological sectors. Results from fluctuating ensembles show no dependence on 
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Figure 6: Effective pion mass as a function of time from one external configuration 
at k = 0.24 
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Figure 7: Effective pion mass as a function of time from ten external configura- 
tions at k — 0.24 
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the topological charge. On the other hand, those projected from non-fluctuating 
ensembles show approximately the same behaviour as simulations completely 
without tunneling. This suggests a rather subtle dynamical effect. 
In Sect. [7| we found that external homogeneous plaquette configurations with 
fixed Polyakov loop values Pt and Pl exhibit completly irregular behaviour. After 
averaging over Pj- and Pl we regain the effective mass results characteristic for 
the topological charge sectors these configurations lie in. 

We conclude that there is a need to obtain a better understanding of the interplay 
between the dynamical mass generation and topological sectors. We would like 
to point out that higher statistics runs could reveal similar phenomena in other 
models with non-trivial topological structure. 
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